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' Abstract. Analysis of both pulse timing and pulse shape variations of the isolated 
l/^ pulsar PSR B1828-11 shows highly correlated and strong Fourier power at periods 

o ; 



^ , ~ 1000,500,250, and 167 d (Stairs et al. |2000| ). The only description based on a free 

I precession of star's rigid crust coupled to the magnetic dipole torque, explains the 500- 

O I component, as the fundamental Fourier frequency, with its harmonic 250-component 
I 



O (Link & Epstein 2001). In this paper, we study a time- varying magnetic field model and 
j/j show that if the dipole moment vector rotates with a period nearly equal to the longest 
(assumed fundamental) observed period (~ 1000 d) relative to the star's body axes, the 



> . . 

!^ ' resulting magnetic torque may produce the whole Fourier spectrum consistently. We also 

X ■ 

$H find the second and fourth harmonics at periods ~ 500 and 250 d are dominant for small 
' wobble angle ~ 3° and large field's inclination angle > 89°. 
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1. Introduction 

The monitoring of long-term and periodic variations both in pulse shape and slow-down 
rate of the isolated pulsar PSR B1828-11 shows strong Fourier power at periods of ~ 
1000, 500, 250, and 167 d, with the strongest one at period ~ 500 d (Stairs et al. |2000| ). 



The close relationship between the periodic changes in the beam shape and the spin-down 
rate of the pulsar suggests the possibility of precession of the spin axis in a rotating body. 
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The precession of the spin axis would provide cyclic changes in the inclination angle x 
between the spin and magnetic symmetry axes. The result will be periodic variations 
both in the observed pulse-profile and spin-down rate of the pulsar. 



Recently Jones & Andersson ( 2001 ) and Link & Epstein (2001) studied a freely pre- 
cessing neutron star due to a small deformation of the star from spherical symmetry 
coupled to a torque such as magnetic dipolc moment, gravitational radiation, etc., and ex- 
plained some part of the observed data. Because of the strong periodicity at period ~ 504 



d seen in the data, Jones & Andersson (2001) reasonably suggested that the actual free 
precession period is Pprc = 1009 d. A coupling between the magnetic dipole moment 
and star's spin axis can provide a strong modulation at period Ppi.c/2 ~ 504 d, when 
the magnetic dipole is nearly orthogonal to the star's deformation axis. But their model 
could not explain the strong Fourier component corresponding to a period of ~ 250 d (see 
Stairs et al. 2000| ). The latter component has a significant contribution in the observed 
variations of period residual Ap, its derivative Ap, and pulse shape. For this reason. Link 



& Epstein (2001) assumed that the strongest Fourier component (~ 500 d) represent 
the actual free precession period. They found that for a small deformation parameter of 
e — {I3 — ~ 9 X 10^^, a free precession of the angular momentum axis around the 

symmetry axis of the crust could provide a period at Pprc — 511 d. Here Ii = I2 < I3 are 
the principle moment inertia of the star. Further, they showed that a coupling of nearly 
orthogonal (fixed to the body of the star) magnetic dipole moment to the spin axis would 
provide the observed harmonic at period ~ 250 d. Their model has good agreement with 
observations in the pulse period, but as they mentioned, it failed to explain the Fourier 
component at period ~ 1000 d seen in the data (as well as 167 d). 

The existence of precession in a neutron star is in strong confiict with the superfluid 
models for the neutron star interior structure. These models have successfully explained 
the glitch phenomena (with both pre- and post-glitch behavior) in most neutron stars 
in which the pinned vortices to the star crust become partially unpinned during a glitch 
(Alpar et al. 1984 ). As shown by Shaham ( |l977|) and Sedrakian et al. ( 1999 ), the preces- 
sion should be damped out by the pinned (even imperfect) vortices on a time scale of few 
precession periods. For example, PSR B1828-11 with typical degree of vortex pinning, 
^pinned /-^star ^ 1.4% (indicated by pulsar glitches in stars that frequently glitch), would 
precess for ^ 40 sec, far shorter than the observed periods (Link & Epstein ^001 ). Here 
^star is the total moment inertia of the star, while /pinnod is the portion of star's fluid 
moment inertia that is pinned to the crust. 

The free precession description provides an effective decoupling between the internal 



superfluid and the crust. Recently, Link & Cutler (2002) studied the problem more care- 
fully by considering dynamics of the pinned vortices in a free precessing star under both 
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Magnus, fm, and hydrodynamics (due to the precession), /p, forces. They found that 
the precessional (free) motion itself prevents the vortex pinning process and keeps the 
vortices unpinned in the crust of PSR B1828-11 while precessing, for a force density (per 
unit length) fp ^ 10^^ dyn cm^^. As a result, they found that partially pinned- vortex 
configuration cannot be static. 

The effective core-crust decoupling causes the core and the crust to rotate at different 
rates (Sedrakian et al. 199£).[| The latter would increase core magnetic flux-tube dis- 
placement relative to the crust, and then sustain the magnetic stresses on the solid crust 



in forcing it to break (in platelets) and move as the star rotate (Ruderman 1991a,l:). The 
stresses are strong enough to move the crustal lattice by continual cracking, buckling, 
or plastic flow to relative stresses beyond the lattice yield strength. As a result, because 
of the very high electrical conductivity of the crust (cr ^ 10^^ s^^), the foot points of 
external magnetic field lines move with the conducting plates in which the field is entan- 
gled. Furthermore, since the core magnetic flux tubes are frozen into the core's fluid, the 
precessing crust drags them and then increases core flux-tube displacements. Therefore, 
during precession of the crust such plate motion is unavoidable. In addition, as shown 



by Malkus (1963, 1968) the precessional motion of the star exerts torques to the core 
and/or crust resulting from shearing flow at the thin core-crust boundary region. These 
torques, so-called precessional torque, are able to sustain a turbulent hydromagnetic flow 
in the boundary region, and then increase local magnetic field strength. This would excite 
convective fluid motions in the core-crust boundary, increase magnetic stresses on the 
crust and cause it to break down in platelets. 

In this paper, motivated by the above conjecture, we suggest that the magnetic field 
may vary somewhat with time, relative to the body axes of the star. The question that 
arises now is whether the whole observed Fourier spectrum of PSR B1828-11, can be con- 
sistently generated by a time-varying magnetic field during the course of free precession 
of the star. In other words, under what conditions will the observed cyclical changes in 
the timing data be produced by precession of the star's crust coupled to the magnetic 
dipole torque of a time-varying magnetic field. In section ^ we address this question in 
detail. Following Link & Epstein ( ^001 ) we assume that the star precesses freely around 
the spin axis, but with period Ppro — 1000 d. Then we show that the magnetic torque 
exerted by a dipole moment may produce the other observed harmonics as seen in data, 
if the magnetic dipole vector rotates with a period close to Ppre relative to the star's 
body axes. Section ^ is devoted to further discussion. 



^ Actually Sedrakian et al. (1999) have shown that even for partially pinned vortices the core 
and the crust would rotate at different angular velocities. 
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2. The precession 

Consider a rigid, biaxial rotating star with the principle axes ei, 62, 63 and corresponding 
principle moment of inertia h = I2 ^ h- The star's angular momentum L is misaligned 
to the symmetry axis 63 by a wobble angle 6, ie. L ■ 63 = LcosO. In general, we assume 
that the stellar magnetic field is dipolar and changes with time as 

m = TOosinx/i(t)ei + mo sinx/2(i)e2 + toq cos x/3(i)e3, (1) 
where toq is the average value of |m| over a period 2t: /ujp and x is the angle between m 
and 63. The functions fi{t),f2{t) and fz{t) are arbitrary functions in time and will be 
determined later by using the data. The equations of motion in the corotating frame are 
I ■ — + a; X L = "^('^ X m) X m - ' m)('^ x m), (2) 

where R is the average radius of the star. The first term of the magnetic torque, Tff , is 
due to the far-field radiation and has components both parallel and perpendicular to the 
spin axis. It is responsible for spinning down the star. The second term, Tnf , represents 
the near-field radiation torque and is exactly perpendicular to the spin axis. It has no 
contribution to the energy/ angular momentum transfer from the star. This torque does 
affect the wobble angle and spin rate of a freely precessing star. Following Link & Epstein 



( 2001 ), for small wobble angle 6' ~ 3°, suggested by the observed pulse shape variations 
of PSR B1828-11 over one precession period, and for small oblateness e ~ 10~®, we have 
(woTff)^^ ^ (woTnf)"^ < 60 <C <C 1. Here cjq is the angular frequency of the star, Tnf 
(~ 10"^ yr), and Tg 10* yr) are the corresponding near- and far-field radiation torque 
time scales, respectively. So up to the first order of 9, the magnetic torques in the RHS 
of equation can be neglected. Therefore in this order, the angular velocity vector a; 
precesses freely around the star's symmetry axis as 

u){t) 2± 9uJo cos(ti)pt + I3p)ei + Oloq sin(ci;pt -f /3p)e2 -|- ^063, (3) 
with the precession frequency ujp = ecjo, and a constant phase [3p. For the case of PSR 
B1828-11, observations suggest that e ~ 4.7 x 10"^ (Stairs et al. |00^. Then ojp ~ 



7.29 X 10"* Hz or equivalently Ppro = 27r/wp ~ 997 d. 
2.1. Timing 

The observed timing behavior can be understood by considering the contribution of other 
torque's components in the variations of the spin rate. By multiplying u; in equations (^ 
we have 

Equation (0) shows the torque-induced variations in the spin rate of star. From equation 
(^ it is clear that the u) ■ T term does not depend on the near-field torque. It contributes 
to the spin rate change only through the negligible final term. Using equations (|l|) and 
(0), by calculating u; • Tff and subtracting equation (0) from the secular spin down of the 
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star (in the absence of precession), — (wsin^ x/'''s){H/ Ii)^ one can find the spin rate due 

to the far-field torque as (dropping constant terms) 
Aw 1 h 



Wo 



Tff II 



cos 



Xfl{t) + esin2x ( cos(wpi + I3p)h{t)h{t) + sin(wpi + I3j^f2{t)h{t) 

-— sin2 X ( sin2(wpi + f3p)fl{t) + cos^uopt + (3p)fl{t) 

-2sm{ujpt + Pp)cos{ujpt + Pp)fi{t)f2it) ^ . 
Using Fourier expansion, we expand the functions /i(t), ,f2{t), and f3{t) as follow 
fi{t) X! cos(riWd + n/3d), 

/2(0 = X! ^" sin(nwd + nPd), 

fsit) ^^Cn cos{nuJd + n(3d), (5a) 

where ujd is frequency of the magnetic field's variation and (3d is constant. The coefficients 
an,bn, and c„ will be determined by fitting the data. Equation (|^) shows that the spin 
rate variations depend on both the precession frequency and the variation frequency of 
the dipolc field. For the case oo = co = 1 and a„ = = c„ = for n > 1, equation 
(^ reduces to one obtained by Link & Epstein ( 2001 ) (except by a factor /3//1). 

To obtain the reported spectrum of PSR B1828-11, it is enough to consider n = 
and 1 terms in equation (^ only. The n > 2 terms will produce the higher harmonics 
which will be discussed later. Simply a correct behavior of the spin rate can be found 
by setting ci — 0, cq = 1, and bi ~ —ai. Therefore equation (|^) reduces to (dropping 
constant terms) 

— ~ — ^ sin 2x f ao cos(wpt + (3p) + ai cos[(wp + ujd)t + (3p + Pd] ") 

Wo TfT ii L V / 

9 / 

-- sin^ X y 2aoai cos(wrft + Pd) ~ 2aoai cos[(2wp + ujd)t + 2/3p + (3d\ 
-al cos(2wpi + 2l3d) -a\ cos[2(wp + Ud)t + 2(/3p + /3d)] ) . (6) 
As one expected, the expression for observable variations in period derivative, Ap, will 
be modified as well. The star's residual in p is owing to both torque variation, equation 
(^, and the geometrical effect. The later is due to the orientation of the star's angular 
velocity vector a; relative to the observer. As expected, the torque effects dominate the 
geometrical effects by a factor {P"^^^! 'kP^ts){Ij,I sin^ X — 100 ~ 1000 for the precession 
period ~ 1000 d, and so we neglected it here. Therefore 

-^Aw^-^ 
27r T 



■ sm 



Q 2x ao cos(wpi + /3p) + oi cos[(wp + ujd)t + Pp + Pd 
X ( 2aoai cos(wdt + /3d) - 2aoai cos[(2wp + Wd)i + 2/3p + Pd 
-Oq cos(2wpi + 2Pd) -af cos[2(wp + ujd)t + 2(/3p + /3d)] 



(7) 

where T = (rff/2sin^ x)(/i//3) 2± tage is approximately equal to the characteristic spin- 
down age and Po is the spin period of star. Equation (Q) gives the period derivative 
residual due to far-field torque variations. The observations showed that both 500 and 
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250 Fourier components are dominant and have comparable amplitudes (Stairs et al. 



200C, Link & Epstein 2001). By assuming the precession frequency as LOp ~ 27r/500 d~^, 
one can get the Link & Epstein's results for |ao| S> |ai| and LOd ^ ^p- But by setting 
u)p + LOd — 27r/508 d^^ (chosen by fitting the data), or equivalently 1/Pprc + ^/Pd — 
1/508 d^^, we have Pd/Ppic — (^pre/508 — 1)"^ where Pd and Ppro measured in days. 
Unfortunately the 1000 Fourier component in the data is not as accurate as the other 
components. So for Ppio ~ 1016 d, we have Pd ~ Ppic or Ud — ujp- The value of Pd will 
be larger (lower) than Pprc, if Ppic (> 508 d) is lower (larger) than 1016 d. Therefore 
equation (0) reduces to (for (3p = and f3d = 0) 



cot X flo cos(27rt/1016) + ai cos(27ri/508) 



-- 2aoai cos(67rt/1016) + al cos(27rt/254) 



(8) 

where t measured in day. Here we ignore the 1000 d and 500 d contributions to the 0^ 
order. For |ao| "C |ai|, the 250-component will be comparable to the 500-componcnt if 
we have {aiO/A) tanx > 1, or tanx > 4/(ai0). For a small 9 (ai > 1), one finds that the 
magnetic dipole moment must be nearly orthogonal to the symmetry axis 63. Hence for 
X > 89° and + ~ 27r/500 d~^, the most dominant terms are the second and forth 
harmonics, 500 d and 250 d, that is in good agreement with observed data. Since the 
proposed inclination angle between star's spin axis and magnetic field's symmetry axis is 
nearly right angle, x > 89°, one may consider the rotation of dipole vector as a magnetic 
poles reversal. We will get back to this point later. 

It is interesting to note that equation (||) includes 1000, 500, 333, and 250 d Fourier 
components. Further, by considering n > 2 terms in equation (^, one can find the higher 
harmonics in Ap. These terms were missing in the Link & Epstein's model. In table 1 we 
compare the time-varying magnetic field model with one suggested by Link & Epstein 



(2001), and the observations made by Stairs et al. (200C) 



3. Discussion 

In this paper, motivated by the effective core-crust decoupling during the precessional 



motion (Link & Cutler 2002), we considered the case of time- varying magnetic field of 
the star. Then we endeavored to find out under what condition a time-varying mag- 
netic field is able to provide a consistent explanation for the reported PSR B1828-11 



timing analysis by Stairs et al. (2000). We studied the free precession of spin axis of 
PSR B1828-11 under the magnetic radiation torque caused by an inclined time-varying 
magnetic dipole moment vector m(i), with constant inclination angle x- In general, we 
assumed that the dipole field vector changes with time relative to the star's body axes. 
Then we expanded its components in terms of the Fourier expansion as mi — m. ■ ei ^ 
"^o sinxI]n=o cosinud + n(3d), ma = m • 63 = mo sinxI]n=o ^" sin(7itjrf + n/3rf), and 
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7713 = m-e3 = mo ainxJ2n=Q'^nC0s{nLUd + nPd), where a„, 6„, and c„ were determined by 
fitting the data. Finally, we showed that if LUd/tOp ~ Pprc/SOS — 1, one may consistently 
explain the whole observed spectrmn of the Fourier power analysis of PSR B1828-11. 
We find an acceptable fit to the data with a precession period of 1015 d, a wobble 
angle 6 — 3°. 2, and the inclination angle x = 89° between the magnetic dipole and 
star's symmetry axis. Note that the chosen Fourier expansion coefficients are ao ~ .01, 
—bi = ai = 1 = Co, and a„ = &o — &n — c„_i = for n > 2. 

The time- varying magnetic field model can also explain the observed timing data for 
PSR B1642-03. The analysis of timing data of PSR B1642-G3, collected over a span of 30 
years, exhibit strong Fourier power at periods ~ 5000, 2500, and 1250 days (Shabanova 



et al. 2001). The suggested wobble and magnetic field inclination angles are 6* ~ 0°.8 and 
X ~ 60°, respectively. Similar to PSR B1828-11, the spectra of PSR B1642-03 show wide 
spectral features at periods ~ 2500 d and 1250 d. The pulse shape variations were not 
detected, probably, due to their small amplitudes. Furthermore, it is interesting to note 
that both PSR B1828-11 and PSR B1642-03 exhibit features around the sixth harmonic, 
6ujp, ie. 167 d and 667 d, respectively. As seen from equation (|^), by including the n = 2 
term in the expansion, one can reasonably get these harmonics in Ap. These terms were 
missing in previous studies. 

The necessary equality relation between Pd and Ppj-c requires the magnetic poles at the 
surface of the star move with relative velocity Vjcv ^ "^TrR/Ppic ^ 7 x 10"^ cm respect 
to the body axes. On the other hand, the magnetic poles reverse every Pprc/2 ^ 500 d. 
Because of the very high electrical conductivity of the solid crust {a ^ 10^^ s^^), this 
result would be hardly acceptable.^ Variation of the magnetic dipole vector with time in a 
neutron star (with solid crust and no active convective zone) may be understood through 
the so-called neutron star crustal tectonics scenario which has been proposed originally 
to explain magnetic dipole evolution and resulting observable features in millisecond 



pulsars, low-mass X-ray binaries and radio pulsars (Rudcrman 1991a,b). The solid crustal 
lattice of neutron star is subject to various strong stresses. The pinned superfluid neutron 



^ Though such a short magnetic cycle has not been observed in neutron stars yet, the early 
observations of A- type stars (a- variables), with kilogauss magnetic field strength, showed large 
amplitude, nearly symmetric magnetic reversals in periods ranging from 4 to 9 days, close to 



the periods of the stars (Babcock 1958). Several recent observations from the young rapidly 
rotating stars confirmed the existence of the solar-type magnetic cycle with Po/Pcyc — 10~* 



(Brandenburg et al. 199^ , Kitchatinov et al. 200C ) . Of course these stars presumably have active 
convection zone, for the case PSR B1828-11 with proposed period for the magnetic cycle, we 
have Po/Pcyc — 10~^, which is smaller by 5 orders of magnitude relative to one obtained for the 
young rotating stars. This may agree with the fact that in neutron stars the convective fluid 
motions are hardly excited. 
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star vortex lines exert a strong force on the crustal lattice of nuclei which pin them 



(Anderson & Itoh 1975, Ruderman 1976, Alpar et al. 1984). Further, the evolving core 
magnetic flux tubes which pass through the crust, pull it strongly at the base of the crust 



Srinivasan et al. (199C). In the rapidly rotating weakly magnetized neutron stars such as 
millisecond pulsars and low-mass X-ray binaries, the lattice stresses from pinned vortices 
are dominant, while in the older pulsars such as radio pulsars with strong magnetic 
fields, the magnetic stresses from core flux-tube displacement may become important. 
For PSR B1828-11 and PSR B1642-03 with magnetic field strength B - lO^^ q ^nd 
effectively unpinned superfluid vortices (Link & Cutler 2002| ), the latter case is more 
appropriate. The quantized magnetic flux tubes in a core's type II superconducting proton 
sea terminate at the base of the crust. These flux tubes move in response to changes in 
the positions of neutron star's core superfluid. If the crust were to remain immobile the 
shear stress, S{B), on the base of the highly conducting crust (from core magnetic flux 
tube motion) could grow to reach 

Here Be > 10^^ G is the average magnetic field in each core flux-tube, and B is the 
average magnetic field through the crust. The maximum stress that crust could bear 
before breaking depends on the lattice shear modulus and is calculated by Ruderman 
( |l991a,h| ) as 5max < 10^^ dyn cm^^ (for most stars 5max - 10^^ - lO^"* dyn cm^^). 
If S{B) > S'max, neutron stars with strongly magnetized cores would break their crust 
continually as they rotate. For PSR B1828-11 with B 5x 10^^ G, we have S{B) > S'max- 
Therefore, one would expect a continuous crust breaking and crustal plate motion in 
this pulsar]^ In addition, the precessing crust drags the core magnetic flux tubes which 
are frozen into the core's fluid. This would increase core flux-tube displacements and 
then increase the characteristic velocity of conducting platelets. A typical characteristic 
velocity of a flux-tube array in the stellar core layer just below the core-crust interface 



is given by Ruderman et al. ( |199^ ) as Vc ^ (t^o/10 Hz)(10i2 G/B)~^10-'^ cm For 
PSR B1828-11, Wc 3 X 10^^ cm s^^ which is much smaller than the proposed relative 
velocity for magnetic poles by our calculations, v^cv ^ ^nR/Ppj-c ^ 7 x 10^^ cm s~^. 
We note that in calculation of Vc the effect of precessional motion of the crust was not 
considered. By including the precessional effects, one may expect the value of flux-tube 
velocity Vc to increase significantly. 



^ It is interesting to note that according to the neutron star crustal tectonics scenario the 
magnetic fields in spinning down neutron stars move to achieve a right angle configuration 



relative to the star's spin axis (Ruderman 1991a,b). This is in agreement with our analysis for 
PSR B1128-11 as we found x > 89°. 
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It is worth noting that the required torque variation may as well be that due to 
the internal torque, by the partially pinned vortices during the precessional motion of 
the crust. The internal torques would arise from different coupled components of the 
star which move with different velocities, eg. the mutual friction torque which arise 
from different velocities of vortex lines and superfluid. These torques are able to sustain 
hydromagnetic shear flows and turbulences in the core-crust boundary, excite the fluid 



convection motions, and cause magnetic field variations (Malkus 1963 , 1968 ). Further, 



they affect the motion of the neutron star crust, for example, by tilting away its angular 



velocity vector from alignment with star principle axis (Sedrakian et al. 1999). To find 
a clear picture of dynamics of magnetic field in a precessing neutron star, one has to 
consider the effect of the internal torques. This is currently under investigation (Rezania 



2002) 



Finally, in this paper we showed that a time-varying magnetic field model is able to 
explain consistently the timing analysis of both PSR B1828-11 and PSR B1642-03, if the 
field's symmetry axis rotates with a rate nearly equal to their precession rates, relative to 
the star's body axes. Unfortunately, at this stage, the large speed of the magnetic poles 
at the surface of the star required by this model is difficult to accept. Further studies 
on the evolution and dynamics of magnetic fields in precessing stars (especially the plate 
tectonics model) seem necessary. These will be left for future investigations. 
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Table Captions: 



Table 1: In this table we compare the time- varying magnetic field model for PSR 



B1828-11 with one suggested by Link & Epstein (2001) and the observed data reported 



by Stairs et al. (2000). Pq, Ppro, e = -Po/-fpio, 6*, and x are the star's period, precession 
period, star's oblateness, wobble angle, and field's inclination angle, respectively. No 
observational information is available for the beam pattern. 



Table 1 



Data/Model Pprc (d) e 6 x 

Data ~ 1000, 500, 250, 167 4.7 x 10"'' ~ 3° > 89° 

Time-varying mag. model ~ 1000^, all harmonics 4.7 x 10^^ ~ 3° > 89° 

Link & Epstein's model ~ 500^*, 250'^ 9 x lO^^ -3° > 89° 



^Fundamental 
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